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ABSTRACT 


The  general  characteristics  of  panel  flatter  at  high 
supersonic  Mach  numbers  are  examined  theoretically.  Linear 
plate  theory  and  two-dimensional  first  order  aerodynamics  are 
used.  The  report  attempts  to  clarify  the  important  role  of 
damping,  the  relationship  between  traveling  and  standing  wave 
theories  of  panel  flutter,  and  the  effects  of  edge  conditions. 
The  solution  procedures  and  general  mathematical  behavior  may 
be  of  interest  in  other  stability  problems  characterized  by  the 
appearance  of  complex  eigenvalues. 
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SECTION  I 
INTRODUCTION 


Partel  flutter  is  the  self -excited  oscillation  of  the 
external  skin  of  a  flight  vehicle  when  exposed  to  an  airflow 
on  one  side. 

This  type  of  aeroelastic  instability  has  received  much 
study  during  the  past  fifteen  years,  both  theoretically  and 

i 

experimentally.  The  early  work  of  Sylvester  and  Baker, 

Nelson  and  Cunningham,^  Fung,^  Hedgepeth,^  Movchan,^  Houbolt,^ 
to  mention  a  few  names,  has  been  supplemented  by  much  recent 
work  on  the  subject  (see,  for  example,  Refs.  7  through  14). 

Today,  a  great  quantity  of  literature  on  panel  flutter  exists, 
and  the  problem  is  reasonably  understood,  although  work  still 
remains  to  be  done  to  better  correlate  theory  with  experiment 
for  certain  panel  configurations  and  Mach  numbers.  Fung, 
in  two  excellent  survey  papers  discusses  the  status  of  the  panel 

flutter  problem.  See  also,  Dowell  and  Voss,^  Bohon  and 

X  A-  1 6  17 

Dixon,  Johns,  Kordes,  Tuovila,  and  Guy,  and  Shirk  and 

18 

Cl sen. 

The  present  report  will  review  the  theoretical  character¬ 
istics  of  panel  flutter  at  high  supersonic  Mach  numbers,  and  will 
attempt  to  clarify  some  of  the  loose  ends  in  the  literature 
regarding  the  role  of  damping,  traveling-wave  versus  standing- 
wave  theories,  and  effects  of  edge  conditions.  It  is  hoped 
thereby  to  clearly  present  the  high  Mach  number  panel  flutter 
problem  and  its  ramifications,  some  of  which  may  not  have  been 
apparent  heretofore. 

The  solution  procedures  and  general  mathematical  oehavior 
may  be  of  interest  also  in  other  stability  problems  character¬ 
ized  by  the  appearance  of  complex  eigenvalues. 


SECTION  II 

BASIC  PANEL  FLUTTER  EQUATION  AND  ITS  SOLUTION 


The  basic  panel  flutter  problem  can  be  formulated  by 
considering  a  flat,  rectangular  panel,  simply  supported  on  all 
four  edges,  and  subjected  to  a  supersonic  flow  over  one  side. 
See  Fig.  1.  The  panel  additionally  is  subjected  to  midplanc 
compressive  forces  and  N^  (lbs/in),  rests  on  an  elastic 
foundation  K  (lb/in^),  and  has  a  structural  damping  Gg 
(lbs-sec/in^) .  The  governing  differential  equation  for  this 
situation  is, 


D 


2. 


3lW  \ 


X  X 


(1) 


7>w 

7>t 


where  D  is  the  plate  rigidity,  Ap-.  is  the  aerodynamic  pressure, 
is  the  density  of  the  panel  material,  and  h  is  the  panel 
thickness • 

The  aerodynamic  pressure  for  high  supersonic  Mach  numbers 
(  M  >  1-7  ),  can  be  reasonably  described  by  the  two-dimensional, 

3 

first-order  theory  approximation, 
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i\ 
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l»l  —  X 


t 
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(2) 


This  expression  assumes  the  pressure  on  the  bottom  side  remains 

constant  at  the  free-stream  value  p  oo  »  Note  also  for  large 
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M,  the  above  reduces  to  aerodynamic  piston  theory  since 


2 


s  m"-  r 


M  and  (  l) 


1  . 


Combining  Eqs.  (1)  and  (2),  and  introducing  the  nondimen- 
sional  coordinates  ^  , t\  ,  T  defined  as 


1  - 
( 

X 

CL 

ii 

r  = 

d)0t 

(3) 

will  result 

in  the  bas  ic 

partial  differential 

equation 

for 

panel 

flutter , 

L  -4- 

o  f a\M 
b 

4-  7 

/  o,\  dW 

+  vt)  ^  + 

a7 

■- 

9r 

4-  7 

+-  TT  1*L 
'u  T-1 

+■  TT  W 

+  TT^r,  —  T 
*  9^ 

ttV^( 

7_  *■>?- 

oA  a  W  _ 

0 

(4) 

where  the  following  nondimensional  parameters  have  been  intro¬ 
duced  , 


'X  = 

T*  J  M1'-  1 

dynamic  press, 
par  am. 

(5) 

** 

H 

*■  ^ 

total  damp,  coeff. 

(6) 

7  - 

Ojrf M(MZ-2)^  ?a  Ca  /a\' 

(  (Ml-0V‘J  ?«  c„W 

aero.  damp,  coeffo 

(7) 

7  < 
a  ~ 

7-  i£i 
(70 

effective  struct, 
damp,  coeff. 

t.8) 

7 

aspect  ratio 

(9) 

JL  = 

_Ka* 
t t*D 

foundation  param® 

(10) 

r.  = 

longit.  compress. 

(U) 

TTi  I) 

param. 

r  — 

1  y  - 

Nly  O} 

lateral  compress. 

(12) 

1 

tr 

param. 

In  the  above  expressions,  the  reference  frequency  U)Q  has  been 
chosen  as 


N 


D 


a' 


(13) 


which  physically  represents  the  lowest  natural  frequency  of  a 
two-dimensional  simply-supported  panel  (  0)  with  no 

airflow,  elastic  foundation,  or  midplane  compressive  forces 
present.  Also,  the  total  damping  coefficient  tjy  is  seen  to 
be  simply  the  sum  of  two  parts,  <.  u  aerodynamic  damping 
coefficient  <^A  and  an  effective  structural  damping  coefficient 
.  The  aerodynamic  damping  coefficient  c^A  ,  first  intro¬ 
duced  by  Houbolt,^  ranges  from  0  to  100  and  is  shown  plotted  in 
Fig.  2  for  different  panel  sizes,  altitudes,  and  panel  materials®* 
This  <^A  is  a  more  convenient  parameter  here,  than  the  commonly 
used  mass  density  parameter  .  The  effective  structural 
damping  coefficient  <^-s  is  a  consequence  of  the  assumed  constant 
structural  damping  Qs  which  can  be  expressed  as, 

5s  =  %;  “i  h  (14) 


The  Mach  number  factor  in  braces  is  often  assumed  to  be  its 
aerodynamic  piston  theory  value  of  1. 
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where  =■  =  2  * (critical  damping  ratio)  of  any  mode 

The  form  of  Eq.  (14)  implies  that  for  any  other  mode 


co- 


the  structural  damping  coefficient  <|j  will  be  given  by 
For  typical  panels,  ^  ranges  from  0  to  0.03  approximately. 
The  consequences  of  using  other  values  of  cj.j  for  the  higher 
modes  is  explored  in  Section  V. 


LO\ 


Thus,  the  basic  partial  differential  equation  for  panel 
flutter,  Eq.  (4),  is  seen  to  depend  here  on  essentially  six 
nondimensional  parameters  manely  )\  ,  <frT  ,  Q/j-,  ,  %t.  ,  Fx  ,  ry  .  This 
equation  is  to  be  solved  subject  to  the  simply-supported 
boundary  conditions , 


at  ^  -  0;  1 

.t  \  -  0  ,  1 


*  W  =  O  , 


i*  =  o 


*  V  =  0. 


-  o 

1)7^- 


(15a) 

(15b) 


The  solution  procedure  for  Eq.  (4)  begins  by  seeking 
solutions  in  the  form, 


w(?,-n,r) 


(16) 


where  in  general , 


0  =■  CX  +-  l  tO 


(17) 


The  7\  variation  in  Eq.  (16)  satisfies  the  boundary  conditions 
at  ^  =0,  1.  Unstable  solutions  will  occur  if  d  is  real 
and  positive  (static  instability)  or  complex  with  a  positive  real 


part  (dynamic  insuability).  Inserting  Eq.  (16)  into  Eq .  (A) 
will  yield  the  ordinary  differential  equation. 


where 


(19) 


A  =  A 

(20) 


Br  i-  i.  Br  = 


(21) 


To  obtain  a  good  insight  into  this  problem,  the  basic  ordinary 
differential  equation,  Eq.  (18),  subject  to  the  boundary  con¬ 
ditions,  Eq.  (15a),  will  first  be  solved  thoroughly.  Then, 
using  Eqs.  (19),  (20),  (21),  the  pertinent  physical  parameters 
9S  ;  a/b  2  4.  ;  rx  ,  ry ,  0  will  be  related  to  the  general 

coefficients  C  ,  A  ,  BR  ,  Br  and  physical  panel  flutter 
solutions  will  be  obtained 
- 

The  mode  parameter  m  is  usually  taken  as  rr\  =  1  .  Actually 
all  results  come  out  in  terms  of  wcu/t,  which  can  be  inter¬ 
preted  as  an  effective  aspect  ratio  <x/y'  where  =  b /m  . 
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The  general  solution  of  Eq.  (18)  is, 


_ .  z<  *  |  ^  t 

w(^)  -  c,  £  +  c2  £  ^  c3e  +  c4e  (22) 

where  the  C-m  terms  are  arbitrary  complex  constants,  and  the 
2.™  terms  are  the  four  roots  of  the  complex  characteristic 
equation  of  Eq.  (18).  Inserting  W  into  the  boundary  conditions, 
Eq„  (15a),  yields  the  matrix  equation, 


r  ~\ 

1 

L 

i 

1 

C, 

0 

A 

zi 

< 

Cr 

=  i 

0  " 

e' 

C- 

c5 

0 

,2.  *. 
t,  e 

z^e*1 

T- 

3 

_*•  .  *4- 
2-4  c 

Ct 

0 

l  ) 

(23) 


For  nontrivial  solutions  of  the  constants  Cm,  the  roots  Zy^must 
make  the  determinant  A  equal  to  zero,  that  is, 


l 

1 

1 

1 

A  = 

2. 

A 

** 

p*’ 

eL 

£^5 

e*4 

1  ^  ! 
2.  e 

i 

z,  e 

z,^ 

T_  ? 

For  a  given  value  of  C 

and  of  A 

,  various  vaL 

can  be  selected, 

and  the  four  roots  H-, 

complex  characteristic 

equation 

of  Eq. 

(18) 

-  0 


(24) 


1 


Then  the  complex  function  x  can  be  evaluated,  where  * 

_ _ (25) 

(2,-Zt)(2,-23)C2,-2*)(2,-2J)(Zi-?<X2;.-?4) 

These  values  of  can  be  plotted  in  the  complex 

plane.  The  combination  of  B^and  B>r  which  makes  Bz)  —  i0 

is  a  solution  (eigenvalue)  of  Eqs.  (18)  and  (15a)  for  the  given 
C,  A  combination. 

Equations  (18)  and  (15a)  were  solved  numerically  by  an 
IBM  1620  computer  using  the  above  procedure.  Since  many 
eigenvalues,  B»R  ,  Bx  ,  can  be  found  for  each  C,  A  combination, 
it  was  necessary  to  trace  out  continuously  the  proper  eigenvalue 
branch  by  increasing  A  continuously  from  zero  for  a  fixed  value 
of  C.  It  is  to  be  noted  that  for  A  =  0,  the  eigenvalues  are 
real  (  Bx  =■  o)  ,  and  are  easily  evaluated  analytically  from 
Eq.  (18)  as 


Bk  -  -<17.4-11  n4”  +  8.8647  n  C  (26 

where  n  is  any  integer.  The  eigenvalues  remain  real  until  a 
certain  value  o^  is  reached  above  which  they  become  complex 
(  if  O  j  .  Figures  (3a~3c)  show  the  real  eigenvalues 

(Bl  ~  o)  ,  while  Figs.  (4a-4c)  show  the  complex  eigenvalues. 
Only  the  most  critical  eigenvalue  branches  for  this  problem  are 


'The  function  $  ,  rather  than  the  determinant  A  itself,  is 
evaluated  to  prevent  the  possibility  of  repeated  roots  causing 
the  determinant  to  approach  zero.  Also,  the  function  5?  , 
unlike  the  determinant  A  ,  will  preserve  its  sign  if  one  replaces 
z,  by  ,  etc. 
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indicated  (largest  Br  for  a  given  C,  A  combination) .  Figure  3 
indicates  the  minimum  value  of  A  for  the  first  appearance 
of  complex  eigenvalues. 

Returning  to  the  physical  panel  flutter  problem,  it 
remains  to  relate  the  general  coefficients  C,  A,  BR  ,  £>i  to 
the  pertinent  physical  parameters  A  /  o ^  )  n/b  ;  Jz.  ,  ,  fy  ,  Q 

of  the  problem.  These  relations  are  given  by  Eqs.  (19),  (20), 
and  (21) 0  Since  0  is  to  be  determined  for  a  given  configura- 


tion,  it  is  conveni  ,-nt  to 

rewrite  Eq0  (21)  as 

e  -  <jTe 

— 

(Qr.  -  L  Qt)  =  o 

(  '7) 

where , 

70 

II 

7 T4 

,  fir,  -(tff-Jk. 

(28a) 

Qi  = 

Bx. 

TT4 

(28b) 

Equation  (27)  can  be  solved  for  Q  to  give 


e  =  - 


Jr  ± 

?_ 


N 


(if 


Q-e  +  t  Q.X 


(29) 


This  can  be  further  reduced  to 


0  = 


- 

L 


O0) 


9 


I 


where 


(31a) 


(3lb) 


For  a  given  configuration  defined  by  given  values  of  X  ,  > 

q/b^  £z.}  rx  }  r^,Eqs.  (19)  and  (20)  can  be  used  to  find  C  and  A. 
From  the  appropriate  Figs.  3  and  4,  values  of  13Rand  Bt  are 
found.  Then  QRand  Qt  are  evaluated  from  Eqs.  (28a, b) o  Finally, 
0  =  R  +-  l  u5  is  solved  from  Eqs.  (30)  and  (31a,b). 

The  complete  panel  behavior  is  characterized  by  plotting 
the  R +■  i  co  variation  with  increasing  dynamic  pressure  X  . 
Instability  occurs  when  c<  becomes  positive  (static  type  if 
also  i3  ~  Q  ,  dynamic  type  if  also  £3  X-  0  ) . 

Some  typical  plots  are  shown  in  Fig.  6.  For  the  case  of 
no  damping,  <^T  =  0  ,  instability  does  not  set  in  until  after 
two  undamped  natural  frequencies  have  merged  (hence  the  common 
term  "frequency  coalescence  flutter") .  For  some  damping  present 
>  0  ,  the  instability  sets  in  at  a  somewhat  higher  value 

of  X  •  TSiis  occurs  when  R  >0  >  i«e.,  when 


(32) 


By  routine  algebraic  manipulation,  this  criterion  reduces  to 

IX  >  <Jr  (33) 

-Q* 
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Hence  flatter  (  #>  -  0  )  occurs  when 


Gii 


(34) 


At  this  flutter  condition,  it  can  be  shown  that 
and  thus  the  flutter  frequency  will  be  given  by 


F 


R 


COq 


(35) 


Frequently,  only  the  flutter  condition  is  determined.  Howevei , 
the  violence  of  the  flutcer  can  also  readily  be  obtained  from 
Eq.  (30). 


The  deflection  shape  W(,^,r[,T)  for  any  physical 
situation  is  found  by  considering  the  real  part  of  the  right- 
hand  side  of  Eq.  (16).  The  complex  function  W  ( ^  is  given  by 
Eqc  (22)  where  the  roots  Zw,  are  those  for  the  given  situation, 
and  the  complex  constants  are  found  from  Eq.  (23).  Setting 

=  1  and  solving  the  first  three  equations  of  the  matrix 
equation,  Eq.  (23),  for  (  ; c3  gives 

N  m 


where 

C*  -  1- 


(w- 


(36) 


V 

1 


(3/) 
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N 


Nr.N,  - 

Then,  taking  the  real  part  of  the  right-hand  side  of  Eq.  (16,) 
gives 

w(in1r)  = 

(38) 

This  can  be  plotted  for  various  times  during  one  cycle  (^r  =  m) 
to  give  a  clear  physical  picture  of  the  deflection  shape. 

At  a  flutter  condition,  one  has  c*  = 0  in  the  above  equation. 

Some  deflection  shapes  at  flutter  conditions  ( =■  o  ) 
are  given  in  Figs.  (7a-7£)  for  various  C,A  combinations  and 
the  most  critical  eigenvalue  branch  (largest  Bx  for  this  given 
C,  A  combination),  For  A  -0  ,  the  deflection  shapes  are  simple 
sine  shape  standing-wave  types 

=•  j"  iO  tt  t  co<l  Co  V  (39; 


_  _  _ \ 

J  £  ^  wr  o«.  co  r  —  Wj  as* v lot) 

As  the  coefficient  A  is  increased  for  a  fixed  value  of  C,  the 
deflection  shape  changes  from  a  standing -wave  type  at  low  values 
of  A  where  purely  real  eigenvalues  (Bt  =  o)are  present,  to 
a  traveling -wave  type  at  high  values  of  A  where  complex 
eigenvalues  (  Br  o)  are  present.  Also,  the  deflections  tend 
to  be  concentrated  at  the  rear  end  of  the  panel  for  large 
negative  values  of  C.  In  these  figures,  only  the  first  half 
of  a  period  is  shown  since  the  second  half  is  the  negative  of 
the  first  half. 
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The  solution  procedures  used  here  i or  the  panel  flutter 
problem  represent  an  exact  solution  rather  than  a  modal  solution 
of  Lhe  differential  equation  and  lienee  do  not  possess  con¬ 


vergence  difficulties.  These  procedures  are  analogous 

20 

used  previously  by  Dugund}i  and  Ghareeb  for  solving  a 

r  2 

differential  equation.  See  also  Movchan  ’ 


to  those 
related 


Many  other  physical  problems  are  characterised  by  the 
ordinary  differential  equation  Eq .  (18)  subjected  to  boundary 
conditions  Eq.  (15a).  The  general  solutions  g^ven  by  Figs.  3 
to  5  and  the  mode  shapes  by  Fig.  7  can  b.*  used  to  solve  them 
also. 


'  3 
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SECTION  III 
APPLICATIONS 


The  general  theory  presented  in  Section  II  is  now  applied 
to  various  physical  panel  configurations  in  order  to  show  the 
characteristics  of  panel  flutter  and  the  effects  of  the  six  non- 
dimensional  parameters,  Eqs.  5  to  12.  Essentially,  the  flutter 
condition  (  <*  =  0  )  will  be  examined,  and  plots  of  dynamic 

pressure  parameter  at  flutter,  '\f  ,  versus  total  damping 
coefficient,  ,  will  be  given  for  different  <Vb  ,  and 

Y'y  configurations.  Also,  an  example  of  the  complete  panel 
behavior  will  be  given., 

(a)  Pure  Aspect  Ratio  Effects,  a/ b 

For  this  series  of  panels,  one  considers  Jz  =  r*  —  ry  =•  C5 
i.e.  ,  no  elastic  foundation,  no  longitudinal  compress  i.on  forces 
and  no  lateral  compression  forces  present*  For  these  cases,  only 
djmamic  type  instability  is  possible.,  Figure  8  shows  the  dynamic 
pressure  parameter  at  flutter  'X?  ,  versus  damping  coefficient  o^r 
for  different  aspect  ratios  a/b  •  The  flutter  dynamic  pressure  'Xp 
becomes  large  for  low  aspect  ratios  (high  Q/b  ).  Also,  "Xf 
becomes  independent  of  damping  at  low  values  of  <frr  and 
roughly  proportional  tc  at  high  values  of  .  This  indicates 
a  change  of  panel  flutter  from  a  constant  dynamic  pressure 
phenomenon  at  low  values  of  damping  to  p  constant  velocity 
phenomenon  at  high  values  of  damping  (thin,  light  panels  in  dense 
air).  This  also  permits  one  to  use  the  "static  airforce 
approximation"^,  for  <^r  <1 

The  flutter  frequencies  oJF  are  also  indicated  in  Fig,  8. 
The  deflection  mode  shapes  for  the  points  marked  with  a  heavy 
dot  are  given  by  Figs.  7a,  b,  c  (  °-/b  '  0  l  XF  =  37o;  ic  oo  f  2.0,000') 
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F i gs .  7c,  f  1  a 


( 


4  ;  ,  "  /  ,  and  Fig.  ,'g  1  t 

'Xf  GOtor,o  J .  The  modes  are  seen  to  change  from  slaaoin; 
wave  types  at  low  values  of  damping  <^T  to  traveling  wave 
types  at  high  values  of  .  Also,  the  modes  become  of  very 
short  wave  length,  the  deflections  tend  to  be  concentrated 
at  the  rear  end,  and  the  flutter  frequency  becomes  high  at  large 
values  of  damping  and  low  aspect  raLios  (high  a/b  ), 


(b)  Pure  Elastic.  Foundation  Effect,  k 

For  this  series  or  panels,  one  considers  a/b  -  =  v\,  =  0  1 

i.e. ,  two-dimensional  panels  with  no  longitudinal  or  lateral 
compression  forces  present „  For  these  cases,  again  only  dynamic 
instability  is  possible.  Figure  9  shows  versus  for 
different  elastic  foundation  parameters  k.  The  flatter  dynamic 
pressure  parameter  increases  with  and  with  k„  The 

presence  of  some  little  damping  is  very  important  at  high 
values  of  k,  since  it  can  raise  the  flutter  parameter  9\r 
well  above  the  <^T - o  value  of  ~  343.  Again,  the  flutter 
phenomena  changes  towards  a  constant  velocity  rather  than  a 
constant  dynamic  pressure  phenomenon  as  or  k  becomes  large 

The  flutter  frequencies  ot)F  are  also  indicated  in  Figo  9. 

At  high  values  of  k,  the  %  becomes  the  simple  natural  frequency 
of  the  section  mass-on-spring  foundation.  The  deflection 
mode  shapes  for  the  points  marked  with  a  heavy  dot  are  given  by 
Figs.  7a,  b,  c  (any  value  of  k  ;  Xf  —  370,  2,000,  20,000). 

Again  the  modes  change  from  standing  wave  types  aL  low  values 
of  ce  and  k  to  traveling  wave  Cvpes  at  high  values  of  o  and  ku 
Also,  the  modes  become  of  ver>  short  wave  length,  the  deflections 
tend  to  be  concentrated  at  the  rear  end,  and  the  flutter  frequency 
high  at  large  values  of  and  k^ 


j 


becomes 


(c)  Pure  Longitudinal  Compression  Effect,  r_. 

For  this  series  of  panels,  one  considers  ~  JL.  —  fy  =■  0, 

i.e.,  two-dimensional  panels  with  no  elastic  foundation,  and  no 
lateral  compressive  forces  present.  Figure  10  shows  'Xp  versus 
<^T  for  different  longitudinal  compression  forces  rA  (negative 
rx  indicate  tension  forces) .  The  dynamic  pressure  at  flutter 

i 

'Xf  increases  with  increasing  tension  (negative  fx  )  and  also 
with  increasing  damping  ,  Again  the  flutter  phenomenon  changes 
towards  a  constant  velocity  rather  than  a  constant  dynamic 
pressure  phenomenon  as  <%r  becomes  large.  For  compressive  forces 
rx  >  +  1  ,  static  instability  may  also  occur.  The  nature 

of  these  curves  for  positive  ^  is  better  illustrated  by  a 
cross  plot,  Fig.  11,  which  shows  'X  for  instability  plotted 
versus  rx  .  The  regions  of  dynamic  and  static  instability  are 
readily  apparent  here.  The  point  ^ =  o  } rx = l  represents  the 
Euler  buckling  load  of  the  panel.  It  is  to  be  noted  that 
aerodynamic  forces  may  stabilize  an  otherwise  statically  unstable 
panel. 

The  flutter  frequencies  are  also  indicated  in  Fig.  10. 
The  deflection  mode  shapes  for  the  points  marked  with  a  heavy  dot 
are  given  by  Figs.  7a,  b,  c  (  rx  =  0  ;  XF  -  370  ,  zooo  ,  zo.oooj, 
Figs.  7e,  f  (  rA  -  -3l  ;  =  4 loo,  co.oso^and  Fig.  7j  (rx  =  +5"  ;  'X F  =  37o). 

Again  the  modes  change  from  standing  wave  types  at  low  values 
of  <^T  to  traveling  wave  types  at  high  values  of  ^ .  Also, 
the  modes  become  of  very  short  wave  length,  the  deflections 
tend  to  be  concentrated  at  the  rear  end,  and  the  flutter  fre¬ 
quency  becomes  high  at  large  values  of  and  negative  . 

(d)  Combined  Compression  and  Aspect  Ratio-Effects 

For  this  series  of  panels,  one  considers  o  ,  i.e.,  no 

elastic  foundation  present,  but  with  various  combinations  of  the 
longitudinal  compressive  force  ,  lateral  compressive  force  fy 
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and  aspect  ratio  a/b  .  For  this  categorv,  two  cases  were  con¬ 
sidered  ,  (a)  rx  variable,  ry  -  O  a/b  -  2.  and,  (b)  variable, 

?  ! b  =  2.  . 

For  Case  (a),  the  results  are  shown  in  Figs.  12  and  13. 

The  results  are  similar  in  nature  to  those  given  by  Figs.  10 
and  11,  except  that  now  a  zero  flutter  dynamic  pressure  con¬ 
dition  occurs  at  f\  =  +  l3  which  is  before  the  onset  of  static 
instability  at  rx  =  +  14  .  This  ^Xf.  =  0  condition  can  readily  be 
removed  by  addition  of  a  small  amount  of  damping  <grT  .  In 
fact,  it  can  be  seen  from  Fig.  12  that  this  anomalous  zero 
flutter  dynamic  pressure  condition  merely  implies  flutter  occurs 
at  a  constant  velocity  rather  than  at  a  constant  dynamic  pressure 
here.  The  deflection  mode  shapes  for  the  points  marked  wi -a  a 
heavy  dot  are  given  by  Figs.  7a,  b,  c  (V*  =8  ;  "A P  “  no,  zooo,  2o,ooo)f 

Figs.  7i ,  j  (rx  -  13;  =  xoj  37o)and,Figs.  7k,  &  (r*=  iff  ■  AF  =  2 oo  ,  zooo'). 

The  same  general  remarks  as  for  the  pure  longitudinal  compression 
effect  described  in  Section  IIIc  apply  here  as  well. 

For  Case  (b)  where  ,  the  results  are  shown  in  Figs. 

14  and  15.  The  results  are  similar  in  nature  to  those  given 
by  Figs.  10  and  11.  Static  instability  sets  in  for  +  S' 

The  deflection  mode  shapes  for  the  points  marked  with  a  heavy 
dot  are  given  by  Figs.  7b, c  (  =  S  ;  =  Zooo,  zo, ooo). 

The  same  general  remarks  as  for  the  pure  longitudinal  compression 
effect  described  in  Section  IIIc  apply  here  as  well. 

(e)  Boundary  Support  Effect 

For  this  series  of  panels,  one  considers  Jk.  =  rg<  =  fy  —  0; 
i.e.,  no  elastic  foundation,  no  longitudinal  or  lateral 
compressive  forces  present,  but  with  different  aspect  ratios 
and  different  boundary  support  conditions  on  the  front  and  rear 
edges  of  the  panel.  For  these  cases,  the  basic  Eq .  18  must  be 
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solved  again  subject  to  the  different  boundary  conditions 

present.  Such  calculations  were  performed  for  clamped- 

22 

clamped  panels  by  Movchan  and  for  c lamped-free  panels  by 
Dugundji  and  Ghareeb  .  The  complex  eigenvalues  BR (  Bx  for  the 
clamped-c lamped  case  as  obtained  by  Movchan  are  presented  in 
Fig.  16  along  with  the  simply-supported  ones  used  here  for 
comparison.  The  resulting  plots  of  dynamic  pressure  parameter 
at  flutter  "\F  versus  damping  <^T  are  shown  in  Fig.  17.  In 
general,  the  clamped-clamped  case  gives  a  higher  flutter 
dynamic  pressure  than  the  simply -supported  case.  However,  at 
either  low  aspect  ratios  (high  °/b)  or  at  high  values  of 
damping^,  the  two  AFJs  approach  each  other.  This  is  probably 
due  to  the  shorter  wavelengths  present  in  the  flutter 
deflection  shapes  here,  and  hence  a  lesser  influence  of  the  end 
boundary  conditions  (see  flutter  deflection  shape  discussion  in 
Section  Ilia) .  Also  shown  for  comparison  is  a  clamped-free 
beam  from  results  of  Ref.  20.  Thj s  also  approaches  the  simply- 
supported  case  at  high  values  of  <^r  where  the  mode  shapes  show 
short  wavelength  traveling  waves  present. 

(f)  Complete  Panel  Behavior 

For  interest,  the  series  of  panels  examined  in  Section 
Ilia  (Pure  Aspect  Ratio  Effect)  was  re-investigated  to  give  the 
complete  panel  behavior  instead  of  merely  the  flutter  condition. 

For  this  series  of  panels,  one  has  =  O  =  *~y  =  0  • 

Then  making  use  of  the  theory  of  Section  II,  the  variation  of 
0  -  c(  +  i.  io  versus  the  dynamic  pressure  parameter  A  was 
determined  for  different  aA  and  configurations.  Instead 

of  plotting  the  c*  and  io  versus  ,  a  new  parameter,  the 

amplification  ratio,  A.R. ,  was  introduced.  This  A.R.  is  defined 
as  the  ratio  of  amplitudes  during  one  cycle  of  oscillation,  and 
indicates  the  violence  of  flutter.  It  is  given  by 
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A.R. 


(40) 


■2.TT 


=  e 


c< 

u) 


Figure  18  shows  the  amplification  ratio  A.R.  versus  the 
dynamic  pressure  parameter  ,  for  °/b  = 0}  4-  and  several 

values  of  <^-T  .  Flutter  (  A.R  >  i  )  is  seen  to  set  in  very 
sharply  for  a/b  =  0  and  low  values  of  damping  t  .  For  low  aspect 
ratios  (a/L  ”  and  also  for  high  values  of  c^r  ,  the  flutter 
condition  comes  in  more  mildly.  Also  shown  on  this  plot  are 
some  values  of  the  frequency  £5  associated  with  these 
amplification  ratios. 
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SECTION  IV 

TRAVELING  WAVE  ANALYSIS 


; 


For  low  aspect  ratio  panels,  the  length  along  the  stream 

direction  is  much  larger  than  the  width  (  »  1  )  .  One 

might  consider  such  a  panel  as  an  infinitely  long  strip  of  1 

finite  width  b,  and  seek  traveling  wave  solutions  of  the  basic 

partial  differential  equation.*  Although  the  use  of  first  order 

aerodynamics,  Eq .  2,  has  certain  limitations  when  applied  to 
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traveling  waves  '  ,  it  will  be  used  anyway  in  order  to  assess 
the  differences  between  the  traveling  wave  analysis  and  the 
finite  panel  analysis  of  the  same  mathematical  equation. 

The  governing  differential  equations  for  panel  flutter 
at  high  supersonic  Mach  numbers  is  given  as  before  by  Eqs.  1  and 
2.  Combining  Eqs.  1  and  2,  and  introducing  the  new  nondimen- 
sional  coordinates  }  ,  n 1  ,  T  defined  as 

%  =  y  n  =  ^  r  =  w„t  (4i) 

will  result  in  the  partial  differential  equation  for  panel 
flutter , 

T  3r 
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This  traveling  wave  approach  for  low  aspect  ratio  panels  was 
investigated  by  Dowell  (Ref.  23)  using  the  complete  linearized 
aerodynamic  theory.  Also,  this  traveling  wave  approach  is 
often  used  in  problems  of  cylindrical  shell  flutter  (Ref.  10) 
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where  the  following  new  nondiraensional  parameters  have  been 
introduced , 
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The  above  nondimensionalization ,  Eq.  41  and  parameters  Eqs.  43- 
49  are  now  based  solely  on  the  width  b  for  this  infinite  traveling 

‘•v 

wave  analysis.  The  reference  frequency  U)0  is  here  chosen  as 


0), 


=  U 


D 


(50) 
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The  basic  partial  differential  equation,  Eq .  42,  is  equivalent 
to  the  previous  equation,  Eq .  4. 


The  traveling  wave  solutions  of  Eq .  42  are  sought  in  the 


form 


If  (-t-x) 


-  w0  e 


(51) 


where  %  is  the  wavelength  and  c  is  the  wavespeed  which  in  general 
may  be  complex,  i.e., 


O  —  CR  +■  L  Ox  (52) 

Generally,  it  is  assumed  that  unstable  solutions  will  occur  if 
the  wavespeed  c  has  a  negative  imaginary  part.  The  v\,  variation 
in  Eq.  51  satisfies  the  boundary  conditions  at  —  o ,  1 
Inserting  Eq.  51  into  Eq.  42  will  yield  the  algebraic  equation. 


(53) 


where 


(54) 
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(55) 
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When  no  air  forces  and  no  structural  damping  are  present., 
y  -  o^r  =  0  ,  the  above  expression,  Eq.  53  gives, 

C  =  ±  S  00  (56) 

For  a  panel  with  no  spring  foundation  and  no  midplane  compressive 
forces,  =  0,  the  minimum  value  of  occurs  for 

rr>  =  1  and  £/zb  =  1  and  is  given  as  5  •  ~  1  .  Hence, 

the  reference  wavespeed  co  can  be  interpreted  physically  above 
as  the  vacuum  wavespeed  for  traveling  waves  with  vn=i,  whose 
half  wavelength, Jl/z  ,  equals  the  width,  b  .  This  also  happens 
to  be  the  minimum  vacuum  wavespeed  possible  for  such  a  panel. 

Note  also  the  simple  relation  of  C0  to  the  speed  of  sound  in 
the  material,  CM  ,  Eq.  55. 

Solving  Eq.  53  for  the  wavespeed  c.  in  the  presence  of  air 
forces  and  damping  gives 


c  =  l  111  H 
Co  Sb 

This  can  be  further  reduced  to 

I 

l  -  [Mr*i 
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where 


&.{r] 

cJUjrj 


•x  £  l 

It  n3  b  ^[p] 


(58a) 


(58b) 


The  complete  panel  behavior  is  characterized  by  plotting  the 

A/ 

Gr  f  i  Ci  variation  with  increasing  dynamic  pressure  for 
various  wavelengths  Z/z b  •  Instability  is  assumed  to  occur 
when  c1  becomes  negative.  This  occurs  when 


By  routine  algebraic  manipulation,  this  criterion  reduces  to 


'A  > 


(61) 


Hence,  flutter  (  -  0 j  occurs  when 


/\  =  ZTT3  5 


or  equivalently, 


(62) 


(63) 
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as 


where  $  is  a  function  mainly  of  wavelength  ratio  ^/ab 
defined  by  Eq.  54.  At  the  flutter  condition,  Eq.  59b  yields 
UU\  [f]  =  ±  ^  / 8  k  and  tne  corresponding  wavespeed 

and  frequency  at:  flutter  become, 


(64) 


CO 


F 


(65) 


For  a  panel  with  no  spring  foundation  ar.d  no  midplane  compressive 
forces  (  A.  =  rx  --  fy  ~  0)  ,  the  minimum  value  of  p  is  ~ 

and  occurs  for  m  =  1  and  =  1.  Placing  these  values  of  5 

and  fc/zb  into  the  above  expressions,  Eqs.  62-65  gives  the 
condition  for  the  first  onset  of  flutter,  for  such  panels. 

Figure  19  shows  a  plot  of  the  dynamic  pressure  parameter 

r  sS 

at  flutter  \F  versus  the  damping  <^T  for  the  infinite  panel 
traveling  wave  theory.  Also  shown  for  comparison  are  tne 
results  of  the  previously  obtained  finite  panel  analyses,  Fig.  8, 
when  converted  to  the  parameters  \  f  and  fusing  Eqs.  43  and  44. 
The  traveling  wave  analysis  gives  a  lower  flutter  dynamic 
pressure  than  the  finite  panel  analysis,  particularly  at  low 
values  of  <^r  .  At  the  higher  values  of  ^  ,  the  agreement 
and  trends  are  better  between  the  two  theor’  ■  .  The  flutter 
frequencies  are  also  indicated  in  Fig.  19  and  the 

agreements  arc  fair.  The  deflection  mode  shapes  for  the  infinite 
panel,  traveling  wave  ana  Lysis,  are  simple  sine  shaped  traveling 
'waves  of  wavelength  &  -  2(3  ,  traveling  at  a  wavespeed  c.  -  c0 
and  having  a  frequency  top  =  2.  ;00.  The  corresponding  deflection 
mode;  shapes  of  the  finite  panel  for  the  points  marked  -y  a 
heavy  dot  are  given  in  Figs.  7e.  f  ( a/y  ~  *■'  ;  ~  -°4  6  0  )  -inn 

Figs.  7g,  h  («/t=  10;  ^rB.o36  ,  -BO  )  • 


They  clearly  resemble  traveling  waves  and  arc  qualitatively 
similar  in  wavelength,  wavespeed,  and  frequency  to  the  infinite 
panel,  particularly  at  high  values  of  cj.T  .  These  finite  panel 
though,  show  large  deflection  amplitudes  towards  the  r-  .it  of  th 
panel,  as  compared  with  the  uniform  deflection  amplitudes  of 
the  infinite  panel,  traveling  wave  analysis.  For  convenience, 
a  comparison  of  the  mode  shapes  from  the  two  analyses  is  shown 
in  Fig.  20  for  the  a/L  =  10  }  <jr  =■  .80  case. 

Summarizing,  it  appears  that  an  approximate  idea  of  the 
flutter  speed,  frequency,  wavespeed,  and  wavelength  can  be 
obtained  from  an  infinite  panel,  traveling  wave  analysis  for 
long,  narrow  panels  at  high  values  of  damping  (light,  thin 

panels  in  dense  air).  However,  the  end  effects  still  play 
important  roles  for  panels  of  ‘Vg  =  10  >  and  any  accurate 
estimation  of  the  flutter  characteristics  and  deflection 
shapes  should  be  made  by  finite  panel  analyses. 
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SECTION  V 

EFFECT  OF  ARBITRARY  STRUCTURAL  DAMPING 


All  of  the  preceding  analyses  have  assumed  that  the  actual 
structural  damping  coefficient  <%]  ,  of  any  mode  (O j  was  related 
to  that  of  the  fundamental  mode  £Oj  by  the  relation, 

_ 

0]  ~  ^  60;  (66) 

Other  relationships  may  be  assumed  or  measured  experimentally. 

For  example,  the  alternate  relationship, 

<57> 

is  commonly  employed  in  standard  V-g  flutter  analyses  in  industry. 
It  is  of  interest  to  examine  the  effect  and  sensitivity  of  these 
other  variations  on  the  previous  results. 

To  study  these  effects,  it  is  convenient  to  solve  the 
panel  flutter  problem  by  modal  methods.  The  basic  partial 
diffential  equation  of  panel  flutter  is  given  by  Eq.  (4).  Modal 
solutions  of  Eq.  (4)  are  sought  in  the  form, 

N 

W(f,v  ,r)  =  £  <^n(r)  mTT7|  (68) 

\  th 

where  ^,Jr)  i-s  bbe  generalized  coordinate  of  the  n  mode,  and 

the  mode  shape,  sin  virr  %  sin  mir  ^  ,  satisfies  the  simply 

supported  boundary  conditions  Eqs.  (na,b)  on  ail  edges. 

Inserting  Ee.  foe)  into  Eq.  (a),  and  applying  GalcrEin’s  meinod 

will  result  after  some  algebra,  in  the  set  of  ordinary  differential 

equations , 
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(69) 


rt  +  S  —  odd 

\  V  4sh 

TT  *  [_j  S1- 


0 


where , 


nj 5  =  i,v 


N 


E 


h 


4 

Y) 


f 


+  A 


+ 


(70) 


The  summation  above  is  taken  over  all  the  s  terms  for  which 
n  +  s  is  an  odd  integer. 

Consider  now,  for  simplicity,  a  two  mode  analysis.  The 
above  set  of  equations  become, 


A 

dr 


+  e 

u  dr 


12L  q. 

Z  TT^ 


=  0 


^  +  f  <*.  *!»-  +  E,  +  8  \ , 

dr2-  0  Jr  #  3  -tr4'  ^ 


=  0 


(71a) 

(71b) 


In  the  second  equation,  an  arbitrary  factor  f  was  introduced 
to  permit  the  possibility  of  changing  the  amount  of  total 
damping  <^T  of  the  second  mode.  To  investigate  stability,  one 
sets 

'JhC'f')  =  6  (?2) 
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i’lie  character i Stic  determinant  of  the  anove  equations  then 
becomes , 


0  +  <JT0  +  U, 

8~X 

3tt4 


- 

3u  + 


6  +■  -f  ^.T  +  E  2. 


=  0 


(73) 


Upon  expansion,  this  gives  the  characteristic  equation, 


—  Z. 


0  +  H3  0  +  0 


h,  e 


h6  =  o 


(74) 


where , 


H  3 
Ht 

H, 


=  ^T(i  +  -f) 

=  t|  *  El  +  i  <J-T 


F  F  4-  l  d'X 

1 —  I  1 —  2- 


3tT 


(75) 


The  roots  9~c(-rlu)  of  Eq.  (74)  are  examined  as  A  increases 
from  zero  for  any  fixed  configuration.  This  gives  the  complete 

Pa  «  n  1  t.  a  l*  A  •  •  '  A  T  A  t»  *  1  H  *  »  /N  A  A  <  «  >■  O  *  »l«»  A  V\  /%  /•»  A  O  /N  A  ■«  1*  ■!  •  »  A 

ai  l  IT  J_  UCUaVlUL  •  J_UO  t-ClL/J-J-JLl-y  UUCUL  O  WUCU  t/ctumco  pWO-Lt-XVC:  » 

To  investigate  only  the  flutter  condition  c L  =  G  ,  one 
sets  0  =  t  d)F  in  Eq.  (74),.  This  equation  then  splits  into 
a  real  and  an  imaginary  part,  namely, 
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_  4- 
60 


F  ~  ft,  -  f<£]  ^  +  [E'E*-  +  d^tT]  =  0  (76a) 


"‘Jrt1’'  f)  UF  *  f  E0  Wf  =  0 


(76b) 


Solving  Eq .  (76b)  gives  the  flutter  frequency, 


WF  = 


60  p 

60,  “ 


AffE, 
\j  i  +  f 


(77) 


Placing  this  into  Eq.  (76a)  gives,  after  much  algebra,  the  dynamic 
pressure  parameter  at  flutter  9<F  as 


X 


8  .'Li 


'ML) 

■  L  +  f  J\j 


(Ex-E.)  * 


(78) 


Equation  (78)  permits  one  to  investigate  the  effect  of  unequal 
total  damping  ( f  1)  on  the  flutter  characteristics.  A  similar 

9 

such  equation  was  presented  by  Bolotin. 


Returning  to  the  damping  characteristics,  it  is  recalled 
that  the  total  damping  is  the  sum  of  an  aerodynamic  damping  c|A 
and  an  effective  structural  damping  .  It  is  convenient  now 
to  differentiate  between  the  damping  in  each  mode,  namely 


It,  = 

(79) 

'Jti  =  +  <3-sx  =  f  <JT 
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whore  Lhe  effective  structural  damp  in;'  roe  ’"1  ic  Lent  of  the 

the  i mode  is  given  from  Eq.  (8)  as 

fa  =  *  %  <8°> 

<JJ  3 

The  aerodynamic  damping  coefficient  <^A  remains  the  same  for  all 
modes  and  hence  there  is  no  need  to  differentiate  here  between 
t> e  modes  in  Eq.  (79).  Using  Eqs.  (79)  and  (80),  one  may  express 
the  factor  f  as 


—  i  +  (>yV) 
%T!  L  +  ^ 


(81) 


Thus,  f  is  seen  to  depend  on  two  nondimensional  ratios,  namely, 


4l 

%A  U)0 


(82) 


\j/-  — 

'  CO, 


(83) 


Also,  for  interest,  the  total  damping  and  the  undamped  natural 
frequencies  of  this  two  mode  system  can  be  expressed  as, 


(84) 


(0o 


(85) 
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For  any  combination  of  ,  fyi.  the  ratios  ^'/^a  and  ^ 

may  be  evaluated  from  Eqs.  (82),  (83),  and  (85).  The  resulting 
values  of  f  and  gT  from  Eqs.  (81)  and  (84)  may  then  be  inserted 
into  the  two  mode  flutter  formula  Eq.  (78)  to  obtain 


Figure  21  shows  a  plot  of  the  factor  f  versus  %si/<^A  for 
various  values  of  ^  .  Also  shown  on  Fig.  21  is  a  plot  of  the 
important  parameter  (  -2.  JT/l  +-  f  )  versus  f.  For  a  given  ''j/>  ,  as 
fs,/^  increases  from  zero,  the  value  of  f  varies  from  f  =  1 
to  the  asymptotic  value  -?  —  f/V*  The  corresponding  value  of 
( Z XT'/ i *-  £  )  decreases  monotonically  from  unity  to  some  other 
asymptotic  value.  Placing  these  results  into  the  formula 
Eq.  (78),  one  can  see  that  because  of  the  factor  (2JT/l*f)> 
the  addition  of  actual  structural  damping  may  actually 

destabilize  the  system,  particularly  for  systems  where  the  Aero¬ 
dynamic  damping  is  small.  The  maximum  amount  of  this 
destabilization  possible  depends  solely  on  the  effective  struc¬ 
tural  damping  ratio  Y*  and  is  given  from  the  asymptotic  values 
of  f  as , 


Af  (wiH' 


ST  y  u  ct*. 


struct . 


(86) 


In  the  case  of  equal  effective  structural  damping  coefficients, 

Y-  1  ,  the  system  is  always  stabilized  by  the  addition  of 

k 

actual  structural  damping  a.  . 


This  destabilization  occurring  upon  the  addition  of. damping 
has  been  pointed  out  by  Ziegler271,  Bolotin  ,  Johns and 
others. 
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The  previous  theory  was  applied  to  a  panel  with  d/g  =  Jk,— 

-  ^  i.e.,  a  two  dimensional  panel  with  no  elastic 

foundation,  and  no  longitudinal  or  lateral  compressive  forces 

present.  For  this  case,  C0L /uJ*  =  1  ,  and  Uz/u),  =  4  .  Two 

types  of  structural  damping  relationships  were  considered, 

namely  a)  ^  for  which  ^  =  4  ,  and  b)  for 

which  y*  -  1  (equal  effective  structural  damping  coefficients). 

Figure  22  shows  the  dynamic  pressure  parameter  at  flutter  Tv  F 

plotted  versus  the  actual  structural  damping  <^t  present  in  the 

panel ,  for  these  two  cases  at  different  values  of  aerodynamic 

damping.  It  is  seen  that  at  = A  ,  the  addition  of  actual 

structural  damping  4  =  .of  will  reduce  the  flutter  dynamic 

pressure  parameter  from  =  2.74  to  =  2isr8  for  the  ^Ti  ~  cft\ 

case,  while  there  is  a  slight  increase  for  the  -  \°ti'  case. 

At  ^a  =  *  >  the  destabilization  for  the  ^  =  <%,  case  is  much 

less  while  again  there  is  a  slight  increase  in  the 

case.  These  curves  of  af  versus  clearly  illustrate  the 

typical  "looping  back"  of  the  V-g  curves  of  the  standard 

flutter  analysis  used  in  industry.  This  "looping  hack"  is 

thus  seen  to  be  a  result  of  unetiual  effective  structural  damping 

coefficients,  and  not  something  associated  with  the  use  of 

n  * 

viscous  type  damping  rather  than  type  damping. 

To  conclude  this  discussion,  it  might  be  interesting  to 
give  a  simplified  expression  for  the  flutter  parameter  TvF  based 
on  the  two  mode  approximation,  Eq.  (78).  If  and  o^T  <<  L 

Eq.(78)  reduces  simply  to 


IL  is  to  be  noted  that  for  no  aerodynamic  or  structural  damping 
(g^.  =  0)  ,  the  dynamic  pressure  parameter  at  flutter  =  2.74- 
for  these  two  mode  analyses  rather  than  the  exact  value  of 
^,-  =  343.  A  four  mode  analysis  should  actually  be  done  for 
numerical  accuracy.  Figure  22  however  does  give  the  proper 
trends . 
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(87) 


A 


F 


This  expression  shows  Che  direct  relationship  of  the  flutter 

parameter  /\p  on  the  factor  (  ^  f)  ,  and  was  first  given 

9 

by  Bolotin.  However,  the  crucial  role  of  the  effective 
structural  damping  ratio  in  establishing  the  correct  value 

of  f  to  be  used  in  this  expression  must  clearly  be  understood 
(see  Fig.  21) . 
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SECTION  VI 
CONCLUSIONS 


This  report:  has  reviewed  the  theoretical  characteristics 
of  panel  flatter  at  high  supersonic  Mach  numbers.  By  using 
linear  plate  theory  and  two-dimensional ,  first  order  aero¬ 
dynamics,  the  problem  was  shown  to  depend  on  six  nondimensional 
parameters,  'X  ,  <^T  ,  Vb  ,  Jk.  ,  rx  ,  ry  . 

An  exact  solution  of  the  resulting  partial  differential 
equations  revealed  the  nature  of  the  eigenvalues  and  their 
general  behavior.  From  these  eigenvalues,  both  static  and 
dynamic  instabilities  can  be  physically  obtained.  The  solution 
procedure  eliminates  difficulties  associated  with  convergence 
of  modal  methods  and  may  also  be  of  interest  in  other  similar 
type  stability  problems,  for  example,  wing  flutter,  flowing 
pipe  lines,  beam  buckling,  etc. 


For  understanding  panel  flutter,  the  effect  of  damping 
is  important.  At  low  values  of  ,  panel  flutter  occurs  at 
constant  dynamic  pressure  q,  and  has  the  appearance  of  standing 
waves.  At  high  values  of  ^-T  (light,  thin  panels  in  dense  air), 
panel  flutter  occurs  at  constant  velocity  V,  and  has  the  appear¬ 
ance  of  traveling  waves.  Many  analyses  have  been  done  previously 
assuming  <^r =  0  (the  "static  air  force  approximation").  This 
is  adequate  in  some  ranges,  but  inadequate  in  others,  particularly 
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implies  flutter  occurs 


at  constant  velocity  V,  rather  than  at  constant  dynamic  pressure 

< 

q  and  does  not  mean  the  system  is  unstable  for  any  airspeed. 


For  pure  aspect  ratio  and  pure  elastic  foundation,  only 
dynamic  type  instabilities  are  possible,  but  with  compressive 
forces  rx  ,  fy  present,  static  type  instabilities  (buckling)  can 
also  occur. 
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The  effects  of  front  and  rear  edge  conditions  on  the  plate 
tend  to  become  unimportant  for  low  aspect  ratios  and  also  for 
high  aerodynamic  damping  <^T,  where  the  resulting  mode  shapes 
begin  to  appear  like  traveling  waves. 

The  flutter  condition  appears  to  set  in  sharply  for 
two-dimensional  panels  at  low  values  of  damping  .  For  low 

aspect  ratios  and  for  high  values  of  <^T  ,  the  flutter  condition 
comes  in  more  mildly. 

Infinite  panel,  traveling  wave  analysis  can  be  used  to 
obtain  an  approximate  idea  of  the  flutter  characteristics  of  low 
aspect  ratio  panels  at  high  values  of  damping  .  The  travel¬ 
ing  wave  analysis  gives  lower  flutter  speeds  than  the  finite 
panel  analysis  and  the  details  of  the  deflection  mode  shapes 
are  somewhat  different,  but  the  general  trends  are  similar 

The  effect  of  adding  structural  damping  to  a  finite 

panel  is  frequently  destabilising.  The  amount  of  this  destabili¬ 
zation  depends  on  the  relative  amount  of  structural  damping 
added  to  each  mode  of  the  panel,  and  is  characterized  primarily 
by  the  ratio  'f'  .  Limits  on  the  maximum  possible  destabilization 
are  established  for  a  two  mode  analysis.  If  the  effective 
structural  damping  is  added  equally  ( ~f~ =  1  ) ,  the  system  is 
always  stabilized. 
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FIG.  I  BASIC  PANEL  CONFIGURATION 
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FIG.  3a  REAL  EIGENVALUES  (Bx  =  0  CASE) 
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FIG.  3b  REAL  EIGENVALUES  (Bx  =  0  CASE) 
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W(f>7)tr)  =  W  (£)  sin  m^rTyjS^*  ia;^T 
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FIG.  6  TYPICAL  PLOTS  OF  PANEL  BEHAVIOR 
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FIG.  7k,  JL  MODE  SHAPES 
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FIG.  14  COMBINED  COMPRESSION  AND  ASPECT  RATIO  (b) 
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FIG.  16  COMPLEX  EIGENVALUES  FOR  OTHER 
BOUNDARY  CONDITIONS 
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FIG.  18  AMPLIFICATION  RATIO  VS. 

DYNAMIC  PRESSURE  PARAMETER 
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FIG.  19  COMPARISON  OF  TRAVELING  WAVE 
AND  FINITE  PANEL  ANALYSES 
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